William I. Fine Theoretical Physics Institute 
University of Minnesota 



FTPI-MINN-07/29 
UMN-TH-2619/07 
September 2007 



Pionic transitions from X(3872) to \cJ 

S. Dubynskiy 

School of Physics and Astronomy, University of Minnesota, 
Minneapolis, MN 55455 
and 

M.B. Voloshin 

William I. Fine Theoretical Physics Institute, University of Minnesota, 

Minneapolis, MN 55455 
and 

Institute of Theoretical and Experimental Physics, Moscow, 117218 

Abstract 

We consider transitions from the resonance X(3872) to the Xcj states of charmo- 
nium with emission of one or two pions as a means of studying the structure of the X 
resonance. We find that the relative rates for these transitions to the final states with 
different J significantly depend on whether the initial state is a pure charmonium state 
or a four-quark/molecular state. 



1 Introduction 



The narrow resonance A(3872) [JJ [2j [31 H] is an extremely interesting hadronic object: its 
mass is tantalizingly close to the D°D*° threshold, w = M(D°D*°) - M(X) = 0.6 ± 
0.6 Me V [5], and the isotopic symmetry is unusually strongly violated in its decays as fol- 
lows from the co-existence of the decay A (3872) — > 7i + 7i^J/ip on one hand and the decays 
A(3872) -> n + n-TT°J/^ and A(3872) -> jj/ip PEj on the other. The suggested interpreta- 
tions of this resonance include a dominantly molecular state [8j [9], [101 [11] of charmed mesons 
(D° D*° + D*° D°) , possibly of the type discussed in the literature long time ago [121 H3], with 
an admixture of pure charmonium [HJ US], a dominantly 2 3 P± state of charmonium [T6l [T7] . 
and a virtual threshold state of DD* [HI QH [201 EQ . 

A quantitative understanding of the internal structure of this peculiar state is a challeng- 
ing task for experimental studies. As has been previously discussed in the literature [UJ [22j 
[16] a measurement of the rates and of the photon spectrum in the decays A (3872) — > DDj 
can provide a useful insight into certain features of the structure of A. In particular, it has 
been suggested [TH [T7] that the observed properties of A(3872) can be, to an extent, mim- 
icked by a 2 3 Pi state of charmonium, so that any 'molecular' admixture would be viewed as 
a secondary effect due to the coupling to the DD* states. In this picture the main available 
indicator of a significant isospin violation in A(3872), the approximately equal rate of the 
decays A(3872) — > pJ/ip — ► n + it~ J/ip and A(3872) — ► ujJ/ip — > n + -n~ n Q J / ip is explained by 
the kinematical suppression of the isospin-allowed transition X(3872) — > uJ/ip- 

In this paper we discuss the transitions from A (3872) to the Xcj charmonium states with 
emission of one or two pions, which can be studied in addition to the observed processes 
A(3872) — > 7T + 7T~J/tf} and A(3872) — > 7T + 7T~7T J/ifj, and which may be instrumental in 
further exploration of the A resonance. The discussed transitions may be accessible for 
experimental observation and may hold the clue to understanding the isotopic structure of 
the A(3872) and of the prominence of the four-quark component in its internal dynamics. We 
argue that the characteristics of such transitions are generally completely different between 
the possible charmonium and molecular components of the resonance A (3872). The rate 
of the one-pion transition relative to the process with two pions is sensitive to the 1 = 1 
four-quark component of the A(3872), while the isoscalar four-quark component should give 
rise to relative rates of two-pion transitions to the Xcj states with different J, which are very 
likely at variance from those expected for transitions between charmonium levels. 

In Sect. 2 we consider the one- and two-pion transitions for the case where A(3872) is 
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treated as a 2 3 Pi state of charmonium. In this consideration we use the standard approach 
based on the multipole expansion in QCD [23] and the chiral properties of soft pionej. We find 
that of the isospin-conserving two-pion transitions by far the strongest should be the decay 
A (3872) — > tttiXci with the transition to Xc2 being heavily suppressed and the transition 
to Xco being forbidden at all in the leading order of the chiral expansion. Given that, the 
rate of the strongest transition in absolute terms is expected to be quite small r(A(3872) — ► 
KKXci) ~ 1 keV. The isospin-breaking transitions A(3872) — > tt°XcJ proceeding due to the 
isotopic symmetry breaking by the u and d quark masses are estimated to be still weaker 
than the two-pion processes. Additionally, the amplitude of the transition A(3872) — > 7r°x c0 
turns out to be zero within the considered approach. 

A significantly different set of one- and two-pion decay rates is expected in the case where 
A (3872) has a substantial four-quark component, so that the emission of the pions can be 
considered as a 'shake off' of light degrees of freedom, as discussed in Sect. 3. In this case no 
strong suppression of a single pion transition should be expected, and in fact such transitions, 
proceeding due to the 1 = 1 part of the wave function of A(3872), are likely to be dominant 
over the kinematically suppressed emission of two pions. Moreover, all three Xcj resonances 
are allowed in the final state, including the Xco, which is additionally favored by the phase 
space factor p 3 . 

At present we cannot reliably predict the transition rates from a four-quark state. How- 
ever a very approximate estimate [H] , based on the general understanding of the hadronic 
parameters, indicates that the single pion transitions A (3872) — > 7i°Xcj should not be very 
strongly suppressed relative to the observed decay A(3872) — ► 7r + 7r~J/7r. Thus one may 
hope that if A(3872) indeed has a significant 1=1 four-quark/molecular component, the 
discussed single-pion transitions may be accessible for experimental studies. 



2 Charmonium 

In this section we treat A (3872) as a charmonium 2 3 P X state, so that the hadronic transitions 
from this resonance can be studied using the standard approach based on the multipole 
expansion. The two-pion transitions arise [23J in the second order in the El interaction with 
the chromoelectric gluon field E a described by the Hamiltonian 

1 A recent discussion of this method in some detail can be found in Ref. [M] , 
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H E1 = - l -enE?{0\ (l) 

where £ a = tf — 1% is the difference of the color generators acting on the quark and antiquark, 
and f is the vector of the relative position of the quark and antiquark. The convention used 
throughout this paper is that the QCD coupling g is included in normalization of the gluon 
field operators. 

The one pion transitions 3 Pi — > 3 Pj tt are induced by the interference of the El in- 
teraction ([1]) and M2 term containing the chromomagnetic field B a and described by the 
Hamiltonian |25j 

H M 2 = eS k Tj (D J B k (0)) a , (2) 

where D is the QCD covariant derivative, uiq is the heavy quark mass, and S = (<j\ + B2) /2 
is the operator of the total spin of the quark-antiquark pair. 

Generally, one also needs to consider the contribution to the one pion transition coming 
from the Hamiltonian of the form: 

H mI = -J&T ^ r i Pi**)" + (DjB k ) a rj L h ] , (3) 
4omQ 

where L is the operator of the orbital momentum of the quark-antiquark pair. However we 
will argue later that this contribution cancels due to the specific form of the gluonic matrix 
element. 

Using the expressions (CQ) and (j2J) one can find the one- and two- pion transition ampli- 
tudes in the form 

AS = A ( 3 P 1 - tttt 3 Pj) = <vnr I E?E%\ 0) A™ , (4) 

A if) = A ( 3 Pi _> n 3 Pj ) = m -l (?r I E a {D . Bk) a + ^ Q) ^ , (5) 

where the heavy quarkonium amplitudes Aff and A^l are defined as follows 

= ±-{l 3 P J \er i gr 3 i a \2 3 P l ) (6) 



32 
and 

42 = T^^ l3pj I £ a r * ^ r ^ + ^ ^ ^ r * ^ 2 3pi > ( ? ) 

with Q being the Green's function of the heavy quark pair in a color octet state. Also in 
the expressions (jSJ) and ((71) we have used the fact that in the matrix elements between color 
singlet states the color factor £ a . . . £ can be replaced by (5 ab /8) £ c . . . £ c . 
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In the leading nonrelativistic limit the spin of the heavy quark pair decouples from the 
coordinate degrees of freedom. Thus, it is convenient to write the amplitude in a form with 
explicitly factorized spin and orbital components. For this purpose we denote Q and r\i the 
spin polarization amplitudes of the initial 3 Pi and the final 3 Pj states. We also introduce the 
polarization amplitudes (pi, ipi and ipij of the initial state and the final state with J = 1 and 
J = 2 correspondingly. The tensor fyj is symmetric and traceless, as appropriate for the 
J = 2 state. All these amplitudes are assumed to be normalized in a standard way, so the 
sums over the polarization states are defined as 

E c;g = E ^ = E fib = E ^ = > 

pol pol pol pol 



E $j ^« = \ (p^ b n + 



pol 

Further we define the initial and the final states 



Su Sjk — - Sij Ski 



3 Pi = e id k Q aj <j) k /V6 , 3 Pj = Pip m bj , (9) 

where a and b stand for the P-wave coordinate parts of the wave functions and the projection 
matrices P-p for the definite total spin J are: 



= 5 tJ /V3 , PW = ey k ^ k /y/2 , PIP = V« • (10) 

Proceeding to the explicit calculation of the heavy quarkonium matrix elements we as- 
sume as in the heavy quark limit that the spin of the heavy quark pair is conserved, i.e. 
(ViCj) = Generally the intermediate states contributing to the Green's function in 
equations (jSJ) and (J7j) could be those corresponding to the S- and D- wave (the P-wave is 
forbidden by the parity), so that two different structures are allowed: (b m | £ a rj Q rj £ a | a n ) = 
S S im 5j n + D (Sij5 mn + S in 5j m — 2/3 5 im 8j n ). However, assuming that the relevant interme- 
diate states are separated by a large energy gap A from the discussed P-wave states, one can 
approximate the Green's function Q as being proportional to the unit operator: Q ~ 1/A. 
In this case the coefficient in front of the S wave is related to that of the D wave: S = 5/3D. 
Thus within the introduced notations the heavy quarkonium amplitudes can be written in 
terms of a scalar quantity A 

*A i&in&jp &ip3jn &ij&np) ^mpr $r i (H) 
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^ijk 2 (,^in$jp ^ip^jn &ij&np) ^qmk^qpr Pyrin 4*r ■ (^^) 

The quantity A depends on details of the heavy quarkonium dynamics and at present is 
highly model dependent, however, it is clear that this quantity cancels in the considered 
ratios of the rates of the single pion and two-pion transitions. 

The gluonic matrix element for the one pion transition in the Eq.(j4]) can be written using 
the form described in [26] 

% (tt I E? (DjB k ) a + (DjB k ) a E%\ 0) = ^ (3 Pj 5 ik - Pi S jk ) , (13) 
with the form factor X related to the well known expression 



X 



TT 



QaQa 



0) = 8vrV2^— '-^Urnl, (14) 
/ m d + m u 



where m u and rrid are the masses of the u and d quarks and f n ~ 130 MeV. It should be 
noted that the specific form of the gluonic matrix element Eq. (fl~3|) . namely the absence of 
the term proportional to 5ij, makes the contribution of the orbital chromomagnetic term ([3]) 
being proportional to the structure (L ■ r) which is obviously equal to zero. 

Using the expressions (jTBl) . (fT4l) and (I9l- (jl~2l) one readily finds the squares of the transition 
amplitudes summed over the polarizations of the final 3 Pj state: 

WT 2 = Pi\A\\ WT 2 = P^\a\\ iWT* = o, (is) 

15 m 2 Q 9m 2 Q 

where p n = \p \ is a pion momentum. We note here that the vanishing of the amplitude 

is a result of our approximation for the Green's function as being proportional to the unit 

operator. The ratio of the amplitudes ([15]) then gives the ratio of the decay rates: 

r 2 : r a : r = 3^ 3 (2) : 5^ 3 (1) :0«1: 2.70 : , (16) 

where p n (i) ~ 334 MeV and p n (2) ~ 285 MeV are the pion momenta in the transitions to 
the 3 Pi and 3 P 2 states correspondingly. 

Let us proceed now to evaluating the two pion transition rates. Using general arguments 
based on the isospin and parity considerations one can see that the pion pair in the transitions 
3 Px — > 3 Pj 7T7T could only be produced in the even partial waves. From the properties of 
the chiral algebra it follows that the expansion of the S- and D-wave amplitudes starts from 
the terms quadratic in the pion momentum (and mass) and therefore kinematically these 
amplitudes are both of the same order in the soft pion limit. 
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Considering the transition 1 ++ — > ++ one can easily see that the pion pair besides the 
D-wave motion in its center of mass frame should also be involved in a D-w&ve motion as a 
whole relatively to the ++ state, thus the discussed amplitude has to be proportional to the 
fourth power of the pion momentum in the chiral expansion. Taking into account the small 
energy release in the transitions 2 3 Pj — > l 3 Pj tttt we should thus expect a strong suppression 
of the transition to the Xco state in comparison with the transitions to the Xci and Xc2- 

The amplitudes of the two pion transitions to the Xci an d Xc2 can readily be evaluated 
using the general form [251 EH] of the amplitude of the dipion creation by the chromoelectric 
field (HTl) 



(tt+tt \E?Ef\Q) = ^(q 2 + ml)6 ij + ^-K[p li p 2j +p lj p 2i -(e 1 e 2 +p 1 -p 2 )S ij } , (17) 

where £1,2 and pi, 2 are the energy and momentum of each pion, q = p\ + p 2 is the total 
4-momentum of the pion pair; b = 9 is the first coefficient in the beta function for QCD with 
three quark flavors, and k is a parameter introduced in [28]. The experimental value of the 
parameter k « 0.2 [29J agrees with the original theoretical estimate [28J. 

The first term in the Eq. (IT7|) arising from the conformal anomaly [21] describes the S 
wave production of the two pions in their cm. frame while the term proportional to k 
corresponds to both S- and P-waves. For an approximate estimate of the amplitude of the 
transition 2 3 Pi — > 1 3 Pi7ttt we neglect the contribution of the n-term since this process is 
dominated by the large S"-wave contribution arising from the first term in Eq. (ll7p . On the 
other hand for the transition 2 3 Pi — > 1 3 P 2 tt7i there is no S'-wave pion pair emission and only 
the .D-wave part of the n-term contributes to the amplitude. In this way one estimates 



^| 2 =150 (^)V + ^) 2 |^! 2 



Aim I ^ — 1 2 I - K 

b 



12tt 



2 



2 



pi pi + - {pi-P2) 2 



2 



^ |2 "f (nr K] ^ 



where in the last transition in the Eq. ffl9l) the averaging on the relative angle between the 
momenta was made. 

For estimates of the rates of 7T7T transitions one needs to evaluate the corresponding phase 
space integrals at the energy release A^ J ) = M(2 3 P l ) - M(l 3 Pj) 



W {J) 
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A numerical integration using the expressions (fIBl) and (|T9|) then yields 

r(2 3p 1 ^ Xc27r7 r)~ iU ' 

so that the transition 2 3 Pi — > % C 2 7T7r is quite small in comparison to the 2 3 Pi — > x c i 7T7t one. 
Clearly, the main factor responsible for such strong dominance of the transition to the Xci 
state is the enhancement by the conformal anomaly of the S wave production of the pions 
in the amplitude described by Eq.( fT7l) . 

Using the expression tHM and the value of the charmed quark mass m c rs 1.4 GeV, one 
can also find the ratio of one- and two-pion decay rates 

r (2 3 Pi -> Y cl vr°) , . 

One can see from this estimate that for the case where X(3872) is considered to be a 
charmonium 2 3 P X state the one-pion transition is significantly suppressed relative to the 
two pion one as a result of small isospin violation. 

The absolute value of the discussed decay rates can be very approximately estimated by 
using the known rate of the decay ip {2S) — > J /ip{lS) n°. Within the similar description [30] 
the gluonic matrix element of the latter decay is found as 
1 5i 

— (Vm 1pm \C r iG r j C Sk | Cn 0n) /3 = A s Sij 6 kmn lj)* m (f) n , (23) 



so that the rate of the decay is given by 

27ml 2vr 



A s \ 2 , (24) 



where the quantity A s is an analog, relevant for the 2 3 Si — > l 3 ^! transition, of the amplitude 
A introduced in Eqs. ljTTj) and (|T2l) . Although A s and A arise from overlap integrals between 
different pairs of charmonium levels, we consider them as being of the same order, A ~ A s , 
in lieu of a more reliable understanding. Proceeding in this way, we estimate 

r (2 3 p -> Xd tt°) ~ \ 3 r W ~+ t°) - °- 06 keV » ( 25 ) 

where p n ~ 574 MeV is the pion momentum in the ip{2S) decay. Even with a large un- 
certainty in this estimate one can conclude that the rate of the strongest of the discussed 
transitions from A(3872) treated as charmonium 2 3 Pi state, 2 3 Pi — > l 3 Pi tttt, is expected 
to be very small. 
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3 Molecular state 



Another possible scenario is to consider the X(3872) as being a four quark state made from 
the two heavy and two light quarks. In this case it appears reasonable to treat the discussed 
transitions X —>■ irXcj and X — > tttiXcJ as a 'shake off' of the light quarks. In particular, the 
spin dependent 'heavy-light' quark interaction is proportional to the inverse power of the 
heavy quark mass rriQ 1 , so that any exchange of the polarization between the light and heavy 
degrees of freedom is expected to be suppressed. Neglecting such exchange the total wave 
function can be written as a product of the two wave functions describing the heavy and 
light degrees of freedom separately. For this purpose we define the polarization amplitudes 

(1 2) (1 2) 

h\ ' and ' associated with the heavy and the light quark pairs in the initial and the 
final states. For the states with the definite J we introduce the polarization amplitudes fa 
and ipi corresponding to the initial and final states with J = 1 and the symmetrical and 
traceless tensor amplitude ipij for the final 3 P2 state. All these amplitudes are assumed to 
be normalized in the standard way and sums over polarizations are the same as in Eq.flS]). 

(1 2) (1 2) 

Finally, introducing the operators L\ ' and H\ ' acting on the polarization amplitudes 
U and hf [hf\H^Hs\hf) ~ 8{j and [lf\LmLn\lf) ~ 5i m 5j n one can write the 
amplitude for the one pion transition 

4 J) = < "Pj \pi e lmn L$ L« HP HP I 3 Pl ) (26) 

with the initial and the final states 

3 A = e ljk hf if fa/VQ , 3 Pj = P^ hf If , (27) 

where the projection matrices P^f are the same as defined in Eq. lflO]) . Using these relations 
one can readily find the amplitudes of the one pion transition to the states with the different 
total spin: 

4°) = 2 Pi fa Ai/V3 , 4 1} = ^ & ^ VhAx/ V2 , Af = -fa Vi Ax ■ (28) 

The quantity A\ accounts for the internal dynamics and at present can only be discussed 
on a model dependent basis. Taking the square of these expressions and summing over the 
polarization amplitudes of the final states complete the routine procedure of finding the 
square of the matrix element 

WT 2 = IpI\m\ iaW 2 =pI\m\ JaW 2 = IpI\m 2 . (29) 
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This relation readily gives us the ratio of the partial decay widths for the one pion transition 

r : r x : T 2 = 4p 3 (0) : 3p 3 n{1) : 5p 3 n{2) » 2.88 : 0.97 : 1 , (30) 

where p v <q\ ~ 436 MeV is the pion momentum in the transition to the Xco state. 

For the two-pion transition from a four-quark state essentially the only guidance is pro- 
vided by general chiral properties, which require that the expansion of the amplitude in 
the pion 4-momenta starts with a bilinear term. In this order only S and D waves for the 
dipion are possible, and one can write the amplitude in terms of the explicit partial wave 
decomposition: 

42 = < 3 Pj | L$ L« ( 6 mn S + $ mn D) I 3 Pl ) j (31) 

where the traceless .D-wave tensor is defined in a standard way: $y = pnP2j +PijP2i — 
2/3 Sij(pi - pi). In the leading order of the chiral expansion the S wave amplitude can 
generally be written as S = a q 2 + bsi£% + cm 2 with a, b and c being coefficients generally of 
order one, so that there is in fact no kinematical suppression of the D-w&ve in comparison 
with the S- wave, and both these terms enter the amplitude on equal footing. It can be also 
noted, once again, that no amplitude for transition to the 1 3 P state arises in this order of 
expansion. 

One finds from the expressions ( l27|) and ( !3~T1) the amplitudes of the two pion transitions 
in the form 

A Xc0 = , A Xcl = (26ij S - $ij D) <f>i i)* A2/V2 , A Xc2 = D e - k fa ip*^i k A 2 , (32) 

where the quantity A2 describes the internal properties of the transition and at present 
is model dependent. Using these expressions and following the standard prescription it is 
straightforward to find the averaged squares of the transition amplitudes 

D| 2 } , (33) 

D\ 2 . (34) 

As previously discussed, in this case there is generally no reason to expect an enhancement 
of the S wave over the D wave, so that the ratio of the two decay rates may be not as 
dramatic as in Eq. fl21l) . 



^ 2 | 2 <j2|5| 2 + ^ 



— *2 — * 2 1 - 1 / - - v 

Pi P2 + - [pi ■ P2) 



A {2) \ 



A2 



P1P2 + - (pi • P2) 
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4 Summary 

The described estimates of the rates of the pionic transitions from X(3872) to the Xcj 
states of charmonium lead us to conclude that these rates exhibit significantly different 
patterns depending on the internal structure of the X(3872). Namely, if this resonance is 
dominantly a charmonium 3 Pi state, the pionic transitions are quite weak with the process 
X(3872) — ► tttiXci being dominant. The single pion transitions, proceeding through the 
isospin breaking by the light quark masses, are even more suppressed with the amplitude of 
the decay X(3872) — > Xco vanishing altogether as a result of a special form of the matrix 
element in Eq. ( 1131) . If however X(3872) is a four-quark/molecular state with a significant 
isovector part in its wave function, the considered single pion transitions should be greatly 
enhanced and exhibit comparable rates for all three Xcj resonances in the final state, as 
described by Eq.(l30"I). One can notice that the yield of the Xco in the final state in this case is 
the highest, while in the pure charmonium case this yield should be strongly suppressed. For 
the transitions with emission of two pions in the case of dominantly molecular X(3872) we 
find that generally there is no reason to expect a very strong dominance of the production of 
the Xci resonance in the final state as compared to x&- The production of the J = state, 
Xcj, is expected to be strongly suppressed in either case due to the general properties of the 
soft-pion expansion. The relative rate of the single and double pion emission is practically 
impossible to estimate quantitatively at present for a molecular X(3872). At a qualitative 
level, once the isospin suppression of the single pion processes is removed by the four-quark 
structure of X(3872), one might expect a larger rate for the one pion transitions than for 
the two-pion processes, given the small energy release in the considered decays. 
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